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Abstract. STAR observes a complex picture of RHIC collisions where correlation effects of 
different origins - initial state geometry, semi-hard scattering, hadronization, as well as final 
state interactions such as quantum intensity interference - coexist. Presenting the measurements 
of flow, mini-jet deformation, modified hadronization, and the Hanbury Brown and Twiss effect, 
we trace the history of the system from the initial to the final state. The resulting picture is 
discussed in the context of identifying the relevant degrees of freedom and the likely equilibration 
mechanism. 



1. Introduction 

This talk is an overview of recent STAR results in the field known as "correlations and 
fluctuations" . Out of multiple reasons to analyze correlations and fluctuations at RHIC known 
before the data became available, there is one that seems to dominate now: by measuring 
correlations and their evolution we learn about equilibration of the system. Is it taking place? 
What is the mechanism? And what is equilibrating? Novel and advanced data analysis 
techniques have been created lately to address these issues. I will focus on the developments 
since Quark Matter 2004. 

Our evidence for equilibration includes elliptic flow, medium modification of mini-jets and of 
charge-dependent correlations. Related issues are the azimuthally-sensitive HBT and the blast 
wave picture. 



2. Elliptic and directed flow as a particle number correlation 

In non-central collisions the initial spatial configuration is anisotropic. But the configuration 
space anisotropy does not produce momentum space anisotropy unless there is another 
mechanism involved, such as local equilibration via rescattering. If the latter is occuring, then 
it is sensible to talk about pressure gradients, and the azimuthal variation thereof. This picture 
naturally calls for analysis of correlations in velocity between local elements of expanding matter, 
but experimentally one is limited to final state particles. The flow analysis in this field began 
with diagonalization of the (momenta-based) sphericity tensor 1 J. In experiments with limited 
acceptance in rapidity at higher energies (AGS, SPS and RHIC) a different technique [2j became 
popular which utilizes the observed azimuthal anisotropy by expanding the particle number 
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Figure 1. Directed flow coefficient as 
a function of pseudorapidity relative to the 
projectile beam in STAR and NA49 



distribution with respect to the reaction plane angle \I/ r into a Fourier series 

d 3 7V 1 d 2 N ~ t \ 1 1 

£ w = s^^ {1+ £ 2t,mC0S[m( *- (1> 

with ^2 typically being the largest component. Cumulant-based analyses, developed later |3j, 
distinguish between two-particle and higher order correlations and separate flow (a multiparticle 
correlation) from the lower-order "non-flow". 

STAR's recent measurements of directed flow (vi), shown in FigEJ use 3-particle cumulants 
to minimize the "non-flow" effects. In the limiting fragmentation picture 4 , what matters for 
the fragmenting system is the relative rapidity with the nearby spectators, hence the choice of 
the abscissa. The data from different energies are consistent with limiting fragmentation. The 
v\ term requires knowledge of \I/ r between and 2tv and thus an azimuthal direction rather 
than orientation of the reaction plane. This direction, given by the spectators, is reconstructed 
event-by-event by means of a position-sensitive zero-degree shower-maximum detector (ZDC 
SMD) which selectively detects fragmentation neutrons behind a RHIC dipole magnet. These 
measurements indicate that the directed flow at RHIC is opposite in azimuthal angle to the 
direction of the spectators (an "anti-flow"). We will return to the velocity aspect of the elliptic 
flow in Sections [SI and [HI 

3. Final state geometry 

While v m characterize the particle distribution in momentum space, observation of the 
corresponding structure in the configuration space rests with the Hanbury Brown-Twiss (HBT) 
interferometry. Herein a correlation function is defined as the ratio of two-particle probability 
density to the uncorrelated reference, which in STAR is obtained by forming mixed pairs. In 
the Cartesian parameterization, 

C fit (q) = 1 + A ex V (-q 2 R 2 - q 2 s R 2 s - q?R 2 ), (2) 

the momenta of a pair are projected onto the transverse plane where one identifies directions 
of their sum and difference. Projections of the momentum difference on those directions 
become, respectively, outward q Q and sideward q s components. The longitudinal component 
qi is orthogonal to both of them. Corresponding to q are the radii parameters R. FigEl shows a 
y/sNN dependence of the quantities representing the best fits of the parameters in EqElto the 
experimental data. The bottom panel presents density of 7r + multiplicity per unit of rapidity, 
which displays a steady logarithmic increase with y/s^N- The latter may be contrasted with 
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Figure 2. Fit parameters for tt~tt~ 
and dN/ dy of 7r + at or near mid- 
rapidity as a function of ^snn. The 
HBT data are either p t integrals or 
represent a p t bin closest to the most 
probable inclusive p t and are taken 
from 7!. The NA49 HBT data are 
preliminary. The dN/ dy data are 
from 0. 



the general lack of energy dependence of the radius parameters which underscores the non- 
triviality of interpreting the HBT results. In the hydrodynamical framework the concept of 
length of homogeneity 1 has been introduced and it is known that the measured radii reflect 
those lengths 2 (or rather, according to 6 j5 a combined effect of the correlation lengths in the 
primordial correlator between the currents emitting the pion radiation and the source size). 
These lengths are smaller than the geometrical source sizes in an expanding (due to pressure 
gradients) and thus inhomogeneous system. A more subtle issue is that of the little relative 
difference between R Q and R S1 which is believed to increase with longer emission. 

The quantitative understanding of this information remains a model-dependent problem. 
In particular, one may assume a local, thermally-equilibrated system and thus ascribe all the 
sources of inhomogeneity (dynamical correlations) to hydrodynamic expansion. Taking this 
approach and fitting the data (including single-particle spectra) with the modified Blast Wave 
model, Retiere and Lisa obtained a set of parameters Remarkably, the azimuthally- varying 
component p2 in their description of the transverse velocity field appears to be non-zero and 
well constrained by the fit. This represents an indirect measurement of the transverse- velocity 
aspect of elliptic flow which complements the particle- number-related V2- 

The azimuthally-dependent radii recently published by STAR (Hj show an effective source 
elongated perpendicularly to the reaction plane. This indicates shorter homogeneity lengths 
(and thus larger gradients) in-plane than out-of-plane, which is consistent with the geometry of 
the early stage of the collision. The early emission is consistent with short emission time. 



Figure 3. Comparison of CD correlations in 77 and <j)\ (a) pp collisions at yfs = 200 GeV, (b) 
peripheral AuAu collisions, (c) central AuAu collisions, (b) and (c) are at ^s^n — 130 GeV. 



4. Charge- dependent number correlations— modified hadronization in the medium 

The charge-dependent (CD) number correlation is defined as a difference between the like-sign 
and the unlike-sign correlations. (In this and subsequent sections we discuss results obtained 
without particle identification). The correlations are defined by comparing same-event pairs 
with mixed-event pairs, obtained from mixing similar events. To characterize and normalize 
correlations, we use 

Ap _ p a ,b ~ PaPb ^ 



where a and b index particles (or in the discretized case, bins in the kinematic space occupied by 
particles), p a ^ is pair density distribution, while p a and p b are the corresponding single-particle 
distributions, whose product forms the reference density p re f. Experimentally one deals with 
fluctuating histogram bin contents n a and n b and EqOl becomes (see jAppendix Aj ) 

Ap n^-n^xn^_Cov[n a ,n b ] 



V n a x rib V n a x n b 

If the variance in the bin population n is dominated by the Poissonian process, Var[n] = n, and if 
in presence of correlations Cov[n a , rib] oc Var[n] (number of correlation sources is proportional to 
multiplicity) this normalization eliminates the trivial multiplicity dependence of the correlation 
amplitude. Eq|H] represents a per-particle measure of correlations. 3 

In paper ^Sj, CD correlations in 77 and (j) have been analyzed simultaneously by using joint 
autocorrelation. FigEl shows a comparison of CD correlations in pp collisions (panel (a)) with 
peripheral and central AuAu collisions. Kinematic cuts on the track separation are applied 
to eliminate the contribution of HBT and Coulomb effects in AuAu collisions (panels (b) and 
(c)) |15j . A much broader HBT correlation for pp collisions contributes the gaussian peak at 
the origin in (a). Transverse momentum conservation also contributes to the structure at the 
origin (the ridge beneath the gaussian peak). The main feature, the negative gaussian peak on 
?7A, is attributed to local charge conservation during hadronization. If the process happens on 
a longitudinally expanding string (as e.g. in Lund model), the strict alternation of positive and 
negative hadrons in rapidity holds (modulo interspersed neutrals). This implies that 77 a for an 
unlike-sign pair tends to be shorter than for the like sign, regardless of 0a, creating a groove 
in the CD plot, clearly seen in Fig^a). The (^-independent groove appears to diminish in 

1 The term "effective length" was used in the original article. 

2 This result was derived for locally thermally equilibrated systems. 

3 In this terminology, the correlation function commonly used in HBT represents a ^er-^ar^c/e-pazr-normalized 
quantity. The choice of normalization depends on the correlation mechanism under study. 



Figure 4. A fantasy event. Positions of spikes indicate 77, (j) 
kinematics of individual particles, their heights - magnitude of 
quantity X assigned to a particle, which could be p t or a derived 
quantity. This event has azimuthal anisotropy, however its V2 may 
be 0, illustrating the difference between number and pt flow. The 
analysis proceeds by selecting a bin of size Srj, 8(f), integrating the 
random field X within the bin, finding the variance of that integral 
in an event sample, then feeding this information into an integral 
equation to yield the correlation of X, as discussed in jAppendix B[ 

amplitude from (a) and (b) to (c). In addition, the negative peak at relative angle changes 
its shape, becoming narrower from peripheral to central AuAu collisions and acquiring more 
rf-(f) symmetry. We hypothesize that this reflects a change in the hadronization geometry from 
one dominated by independent string break-ups in peripheral collisions to a bulk hadronization 
process in which individual strings are no longer relevant. 

5. Charge-independent number correlations — longitudinal minijet deformation 

The charge-independent (CI) number correlation sums the like-sign and the unlike-sign 
correlations. In AuAu collisions, these correlations have been analyzed by STAR in the space 
of pseudorapidity 77 and azimuthal angle (j) [HJ Edj and in the transverse momentum space [121 • 
Our discrete-wavelet technique [IT] extracts information about correlation structure by 
measuring the power spectra of local fluctuations in the density of charged hadrons with respect 
to a mixed-event reference (the so-called "dynamic texture" of the event). As discussed in fT^j . 
this observable is a measure of the gradient of the two-particle correlation function, but takes 
less computing. 4 

The experiment [T^ reveals a reduction in the dynamic texture (and thus in the correlation 
function gradient) along the //-direction at p t > 0.6 GeV/c in the central AuAu collisions, 
compared to an expectation based on peripheral events. HIJING simulations 11 point to 
minijets as the likely source of correlations of this scale and magnitude. The possible mechanism 
of the reduction is a coupling between the minijet fragments and the longitudinally expanding 
bulk medium. In paper [TH], the same correlation structure is investigated in the joint r]-(j) 
autocorrelation technique. It is found that the correlations are broadened in t]a with centrality. 
A reduction in the correlation gradient (" dynamic texture" ) and the elongation of the correlation 
function along 77 a are consistent descriptions of the effect. 

6. pt correlations —minijets and elliptic flow 

We require a direct measurement of correlations in the velocity field, independent of particle 
number correlations. To form an independent measurement, these have to be separated from 
particle number correlations. For a random field (such as particle number or velocity distributed 
in 77 and 0, see Fig0} variance of the integrated content of a bin can be related via an 
integral equation to the correlation function of the field \W\ I17j . The correlation function 
can be reconstructed from the measured variance by solving the integral equation derived 
in \W\ I17j . This process offers a computational advantage since the bin integration involves 
O(N) computations, whereas two-particle correlations require 0(N 2 ) {N is the event multiplicity 
used). By varying the bin size one varies the range of the difference variable in the correlation 
function. 

4 The number of computations required is 0(N), in contrast with 0(N 2 ) required for direct construction of the 
correlation function. 




Figure 5. Top: scale dependence of the reference-subtracted variance A<7p t:n . Bottom: 
reference-subtracted normalized correlation functions reconstructed from A<jp t:n by inversion. 
The percentage intervals indicate centrality for each column-wise pair of plots. 




Figure 6. Same data as in FigEJwith subtracted azimuthal harmonics v\ and V2. 



In order to extract the pt correlation alone, we need to disentangle the effect of n from that 
of pt,i m the variance of p t = Y^ie(rj,^)bmPt,i^ an d from statistical fluctuations. We denote the 
desired variance as Aa^ :n . As an estimator of A(j^ :n , we use Var[p t — np t ]/n — cr? t , where the 
latter term is the inclusive per-hadron p t variance. In the case when p t in the bin, event-by- 
event, derives from the same parent distribution (not the case in FigEJ), A<7p t:n converges to 0, 
despite fluctuations (statistical or not) in n. 

Top panels in FigEl represent scale dependence of the variance Aa^ :n for three centrality 
classes. The bottom panels, representing per-particle normalized correlations plotted on 
difference variables, are obtained from the respective top panels by solving the integral equation 
derived in [HE [El (see jAppendix B| ) which relates variance and correlation. 

While the quadrupole wave in FigEl^bottom row) is the reflection of the elliptic flow in the 
Pt : n correlation function, the peak at angle differences likely reflects minijets |19j . But unlike 



flow, the effect of minijets on the HBT radii at RHIC energy is not part of the standard (local 
equilibrium) picture. Therefore this effect may constitute an under-appreciated factor driving 
the effective radii (FigEJ) down at higher ^snn. 

The amplitude of the minijet peak has a non-monotonic behavior with increasing centrality, 
rising in peripheral and falling in central collisions. The peak is seen to be elongated 
longitudinally in central collisions, showing qualitatively the same behavior as seen with number 
correlations by using wavelet analysis and joint autocorrelation (SectiorJHl)- 

7. Conclusions 

The small effective HBT radii at RHIC indicate that the system is highly inhomogeneous and 
decouples rapidly. STAR has evidence that semi-hard scattering leaves a trace in the "soft" 
Pt domain which is usually considered a safe realm for local equilibrium approximations, thus 
possibly explaining some of the inhomogeneity. Whether or not the minijet correlations present 
a challenge to the picture of local equilibrium which constitutes a basis for much of the HBT 
interpretation infrastructure, will depend on the effective final-state-particle number involved in 
a minijet, or in other words, on the order of those correlations (higher order being more hydro- 
like and thus probably less of a challenge). Therefore a study of higher order correlations and 
cumulants may be a promising future direction. The phenomenon of //-broadening of minijets 
at "soft" p t (Section [3 and can clarify the nature of the medium created in central AuAu 
collisions, provided that the coupling mechanism, apparently sensitive to that nature, is under 
theoretical control. 

Appendix A. Notation 

In our notation i is particle index, n is number of particles within a kinematic cut (bin), (. . .) 
denotes an average over events, p t is an inclusive mean p t per particle, xa (variants: t/a, 0a) is 
difference variable — X{ — 5x is scale (range of local integration, see FigEJ), Ax is the upper 
limit on Sx. Covariance of x and y is denoted as Cov[x,y], variance of x - as Var[x]. 

Appendix B. Prom variance to correlation 

The connection between variance and correlation function - the basis of analysis in SectiorJHl- 
is discussed in detail using integral calculus in [03 and the algebra of discrete bins in [Ej. For 
a random field X(t), the autocorrelation function is defined [IE| as 

p(X, t A ) = X(t)X(t + t A ) (B.l) 

(assuming that correlations do not depend on position t). 

The variance of X, integrated over bin of width 5r], 50, centered at (see Fig^J is 

/ rSrj/2 r 5<f>/2 \ 2 / rdrj/2 M/2 _ \ 2 

Var[X;<fy,50] = / / Xdr)d(j)\ - / / Xdrjdc/)) = (B.2) 

\J-5<n/2 J-S(f)/2 J \J-8rj/2 J-5</)/2 J 
r5ri/2 M/2 r5r)/2 M/2 

/ d Vl / #1 / d m / # 2 [%,?ii)Ife^2)-%^i) x Xfafc)] (B.3) 

J-Srj/2 J-5</)/2 J-Srj/2 J-S^/2 

As we form the reference-subtracted variance Acr^ = Var[X; £77, Sep] — Var[X re f; ^77, 50], the first 
term in the brackets above becomes the reference-subtracted correlation Ap(X, 77A, 0a) and the 
second term cancels with the reference since X = X re f . 

r5ri/2 M/2 rSr]/2 M/2 

Aa 2 x = dru #1 / d m / d0 2 Ap(X,7/i - 772^1 - <fo) (B.4) 

J-5<n/2 J-S(f)/2 J-5r]/2 J -8<j)/2 
rSr] rScf) 

= 2 dr) A 2 d</> A (8r)-riA)(5<t>-<l> A )Ap(X,T) A ,<t> A ) (B.5) 
Jo Jo 



Eg IB. 51 is obtained by making a change of the integration variables from 771,772 to 771 and 
?7A = 771 — 772, integrating over 771, and repeating the same for <j) (see [16 for details). In 
the actual analysis, this integral equation is discretized by partitioning the acceptance into 
ms x n§ microbins of size £ v x e^. X is replaced by (p t — nptj/y/E. This denominator creates 
the reference number correlation p re f (n) in the final expression IB . 61 which makes it a per-particle 
measure: 

ms ' ns Ap(p t : n;ke , les) 

k,i=i yf>ief(n]k£ri,te</>) 
Here the kernel K is a discrete representation of the continuous case : 

(Sri - r] A ) (8<t> - <M) e v £<i>K msns: ki = e v e^(m s - k + ^)(n s - Z + ^) (B.7) 

Knowing the computationally cheaper quantity Acr^ t:n , we solve the integral equation IB. 61 for 
Ap/^/pref using standard numerical techniques |T7| . 
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